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Effects of Rotation on the Revival of 
a Stalled Shock in Supernova Explosions 

Tatsuya Yamasaki^ and Shoichi Yamada^'^ 

ABSTRACT 

In order to infer the effects of rotation on the revival of a stalled shock in su- 
pernova explosions, we investigated steady accretion flows with a standing shock. 
We first obtained a series of solutions for equations describing non-rotating spher- 
ically symmetric flows and confirmed the results of preceding papers that, for a 
given mass accretion rate, there is a critical luminosity of irradiating neutrinos, 
above which there exists no steady solution. Below the critical value, we found 
two branches of solutions; one is stable and the other is unstable against radial 
perturbations. With a simple argument based on the Riemann problem, we can 
identify the critical luminosity as the one, at which the stalled shock revives. We 
also obtained the condition satisfied by the flow velocity for the critical luminos- 
ity, which can be easily applied to the rotational case. If a collapsing star rotates, 
the accretion flow is non-spherical due to centrifugal forces. Flows are accelerated 
near the rotation axis whereas they are decelerated near the equatorial plane. As 
a result, the critical luminosity is lowered, that is, rotation assists the revival of 
a stalled shock. According to our calculations, the critical luminosity is ~ 25% 
lower for the mass accretion rate of lM0/sec and the rotation frequency of 0.1 
Hz at a radius of 1000 km than that of the spherically symmetric flow with the 
same mass accretion rate. We found that the condition of the flow velocity at 
the critical luminosity is first satisfied at the rotation axis. This suggests that 
the shock revival is triggered on the rotation axis and a jet-like explosion ensues. 

Subject headings: stars: rotation — shock waves — hydrodynamics — super- 
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1. Introduction 

Core collapse supernovae play important roles, for example, in star formations, evolu- 
tions of galaxy and the entire universe, and accelerations of cosmic ray particles owing to its 
nucleosynthesis, energetic shock waves and high luminosity. Recently the phenomena are put 
in the spotlight because the gamma-ray bursts are likely to be associated with them. Un- 
fortunately, however, the mechanism of supernova explosion is still unresolved [e.g. Wheeler 
(2003) for reviews]. 

The supernova explosions are commenced by gravitational collapse of massive stars at 
the end of their lives. When the nuclear saturation density is reached at the center of core, 
the bounce occurs, producing a shock wave that starts to propagate outwards in the core. 
It is almost a consensus among the researchers that this shock wave does not yield the 
supernova explosion directly, since the shock loses energy by photo-dissociations of nuclei 
as well as neutrino-cooling and stalls somewhere inside the core. It is widely expected that 
the shock will be revived by the irradiation of neutrinos diffusing out of the proto-neutron 
star, the scenario originally proposed by Wilson (1982). Although many researchers have 
studied this scenario intensively and extensively, the state-of-the-art numerical simulations 
have shown so far that, as long as collapse is spherically symmetric, the stalled shocks do 
not revive (Liebendorfer et al. 2001; Buras et al. 2003; Thompson et al. 2003; Liebendorfer 
et al. 2004). 

These simulations are concerned with the dynamical evolution of the entire core from 
the onset of the collapse till long after the stagnation of the shock. On the contrary. Burrows 
& Goshy (1993) took a different approach, focusing on the revival of the stalled shock wave. 
After the shock is stagnated due to the energy losses mentioned above, the accretion flows 
through the standing shock wave are quasi-steady (see Janka (2001) for the criticism on 
this point). They approximated these flows with steady state solutions for constant mass 
accretion rates and neutrino luminosities. Varying these constant values, they found that 
for a given mass accretion rate there is a critical neutrino luminosity, above which there 
exists no steady solutions. Using this fact, they argued that the revival of stalled shock 
occurs when the neutrino luminosity exceeds this critical value. The existence of the critical 
luminosity was also observed in the numerical experiments done by Janka & Miiller (1996), 
where they artificially increased neutrino luminosities in their failed explosion models and 
saw what happened to the stalled shocks. Merits of these steady state calculations are not 
only that they are computationally simpler but also that they can clearly demonstrate, if 
qualitatively, what is the cause of the failure of explosions. That is, according to this theory, 
an inappropriate combination of the neutrino luminosity and the mass accretion rate. 

So far, we have been talking about the studies of spherically symmetric collapse. How- 
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ever, several observations suggest that the supernova explosions are intrinsically asymmetric 
in general. For example, Leonard et al. (2000) and Wang et al. (2002) observed a few per- 
cent of linear polarizations for photons from collapse-driven supernovae and argued that the 
stellar envelopes are globally asymmetric, oblate or prolate with an aspect ration of ~ 2. 
We also know that supernova remnants are asymmetric in general. We can easily imagine a 
couple of possible causes for this asymmetry, e.g. hydrodynamical instabilities, rotation and 
magnetic fields. In fact, young pulsars are thought to have rather large angular momenta 
(Kaspi et al. 1994; van der Swaluw and Wu 2001) or magnetic fields (Kulkarni 1992), which 
are probably taken over from the pre-coUapse stars. It is also believed that some convective 
motions occur in the supernova core. The large scale non-spherical oscillations of the stand- 
ing accretion shocks may also cause the stellar envelope globally asymmetric(Blondin et al. 
2003). 

Among these effects, we pay particular attention to the stellar rotation in this paper. 
Note, however, that the rotation of the core may not be so rapid after all. Recent evolution 
models by Heger et al. (Hcgcr ct al. 2000, 2003) suggest that the transport of angular 
momentum during the quasi-static evolutionary phase of the progenitor deprives the core 
of substantial fraction of its angular momentum, particularly when the magnetic torque is 
taken into account (Spruit 2002; Heger et al. 2003). If this is really the case, the rotation 
will play no significant role in dynamics of core-collapse as shown by Buras et al. (2003) (see 
also Miiller et al. (2004)). We had better bear in mind, however, that the evolution models 
are based on ID calculations and have some uncertainties in the mechanism and treatment 
of angular momentum transport. 

The dynamics of rotational collapse has been explored numerically by some authors. For 
example, Monchmeyer & Miiller (1989) and Yamada & Sato (1994) demonstrated that rapid 
rotation tends to weaken the prompt shock wave because the centrifugal force prevents the 
core from contracting sufficiently. Fryer & Heger (2000) found that the efficiency of energy 
transfer by convection was reduced by the decrease of effective gravity near the equatorial 
plane due to the centrifugal force and, as a result, the shock did not propagate as much 
as in the non-rotation case. They also asserted that the resultant anisotropy of energy 
transfer would have made explosion jet-like if it had been successful. Kotake et al. (2003) 
estimated the anisotropy of neutrino radiation from the rotational cores. They claimed that 
the deformed neutrino-sphere, oblate in general, yields the neutrino fiux and the heating rate 
that are larger near the rotation axis than near the equatorial plane. Given the results by 
Kotake et al. (2003) and the results by Shimizu et al. (2001) that such anisotropic neutrino 
irradiation could lead to successful jet-like explosion, we may well expect that the rotation 
is helpful for the explosion in the neutrino heating scenario. 
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In order to see the effect of rotation more clearly, if qualitatively, we followed the lead 
by Burrows & Goshy (1993). We solved time-independent Euler equations describing steady 
accretion flows through a standing shock, taking rotation into account. We paid particular 
attention to how the critical luminosity of neutrinos is altered by the stellar rotation and 
demonstrated that the rotation assist the revival of the stalled shock. 

This paper is organized as follows. In the next section, wc discuss the spherically 
symmetric case. We find two branches of steady solutions and discuss the revival of the 
shock wave employing the critical luminosity. We examine in section 3 the stability of 
the two branches of solutions obtained in section 2. In section 4, wc analyze the steady 
axisymmetric solutions with rotation and discuss the effect of rotation on the revival of the 
stalled shock. The final section is devoted to summary and some more discussion. 



In this and next sections, we discuss spherically symmetric fiows, solving one dimensional 
hydrodynamical equations. The calculations in this section are similar to those of Burrows 
& Goshy (1993). After the shock is stalled, mass accretion rates or neutrino luminosities 
change rather slowly and the fiows become quasi-steady. Idealizing this situation, we seek 
in this section steady solutions of the accretion fiow. 

Since our purpose of this paper is not to determine quantitatively the critical luminosity 
of neutrinos, we make following assumptions and simplifications: i) We do not consider 
the fiows outside the shock wave in detail, simply assuming a free fall, and restrict our 
calculations to the inside of the shock surface, ii) Newtonian formulations are adopted 
because the region of our interest is outside of the neutron star and general relativistic 
effects are negligible there, iii) We do not solve neutrino transfer equations assuming that 
the luminosity and energy of the neutrinos are independent of radius. It is also assumed that 
the neutrinos are thermal with a fixed value of temperature (T,^ = 4.5MeV), for simplicity, 
iv) A simplified equation of state and rates of neutrino heating and cooling are employed, v) 
Photo dissociations of nuclei, magnetic fields and convections are neglected, for simplicity, 
although there should exist regions where fiows are convectively unstable. 

Then, the basic equations describing steady spherically symmetric accretion fiows are 
given in the spherical coordinate as. 



2. Spherically Symmetric Flows With a Standing Shock 
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dur ^1 dp ^ GM ^ _ s 

dr p dr ' 



de p dp . 

where Ur, p, p, e denote radial velocity, density, pressure and specific internal energy, re- 
spectively; M, G, M, ^ are mass accretion rate, gravitational constant, neutron star mass 
and net heating rate by neutrino irradiation, respectively As for the neutron star mass, we 
adopt take 1.3M0. 

Pressure and specific internal energy of nucleons, photons, and relativistic particles are 
approximately written as 

Il7r2 . pkT 

_ llTT^ 3 kT 

where k, c, h, mN are Boltzmann constant, speed of light, Planck constant, and nucleon 
mass, respectively; T is temperature of matter and photons. 

As for the heating and cooling processes, we take into account only the absorption and 
emission of neutrinos by nucleons. We adopt the formulae given in Herant et al. (1992) for 
degeneration factor rj = 0; that is 

g = 4.8 X io32 -^'-e(foes s-^) ^, _ ^ ^ ^ 

where Li,^ is the luminosity of electron neutrinos, which is assumed to be equal to that of 
anti-electron neutrinos. 

The outer and inner boundaries of calculations are set at the shock surface and the 
neutrino-sphere, respectively. We impose the Rankine-Hugoniot relations at the shock sur- 
face. The flow into the shock is assumed to be a free fall with the velocity of U{ — —^/2GM/r 
and the temperature of 0. Thus, the outer boundary conditions are written as 



pu^+p^ PfUf, 



(7) 
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where pf is upstream density at the shock surface. The continuity equation is automatically 
satisfied by imposing equation (1) on each side of the shock. We impose the condition that 
the density be 10^^ g/cm^ at the inner boundary; this is approximately corresponding to 
the condition that the optical depth from the neutrino-sphere to infinity equals to 2/3. The 
latter condition was adopted by Burrows & Goshy (1993). 

We note that the position of the inner boundary is determined by the relation 



where a is the Stefan-Boltzmann constant and r^, is the radius of the neutrino-sphere, that 
is, the inner boundary. On the other hand, the position of the outer boundary is determined 
so as to satisfy the above jump conditions. Thus, the problem is regarded as an eigen-value 
problem. We solve these equations for a wide range of values of the mass accretion rate and 
the neutrino luminosity. 

The results of the calculations are shown in figures 1 and 2. In figure 1, we can see there 
are two solutions for a given mass accretion rate and a neutrino luminosity when the neutrino 
luminosity is below a certain critical value. The shock radii of these two solutions differ. 
Prom now on, we refer to the solution with a smaller shock radius as the inner solution and 
to the other solution as the outer solution. As shown in figure 2, as the neutrino luminosity 
is raised with the mass accretion rate fixed, the shock radius for the inner solutions becomes 
larger, whereas that for the outer solutions becomes smaller. Two solutions coincide with 
each other when the luminosity reaches the critical value. For the luminosity over the critical 
value, there is no solution. This critical value of neutrino luminosity is shown in figure 8. 
Note that when the luminosity is close to the critical value, the shock radii are insensitive 
to the luminosity. Although our formulation is different in detail from Burrows & Goshys', 
the results are in good agreement [cf. figure 1 of Burrows & Goshy (1993)]. As mentioned, 
the shock radius for the critical luminosity is almost independent of the mass accretion rate, 
and is about 200 — 300 km according to our calculations. 

This behavior of the solutions is understood if one remind of the adiabatic subsonic 
Bondi-Hoyle fiow (cf. Shapiro and Teukolsky (1983)). Here, we introduce u^osti^)-, as the 
post shock velocity that we would obtain if the shock radius were r. The velocity iipost('") is 
shown as dotted lines in figure 1. The subsonic adiabatic fiows accelerate outside a certain 
radius, whereas they decelerate inside that radius. The critical radius is determined by the 
balance of the gravitational energy and thermal energy and is close to the sonic point in the 
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transonic accretion flow. On the other hand, the free fall velocity decreases monotonically. 
Furthermore, the acceleration of the adiabatic flow is larger than that of the free fall outside 
the critical radius owing to the pressure gradient force. Therefore, if any, there are two radii 
where the accretion flow joins with the free fall stream via the Rankine-Hugoniot relations. In 
other words, the solution curve of the adiabatic accretion flow Ur{r) intersects with 'Upost('") at 
two points. Although the solution curves are somewhat modifled if one takes account of the 
neutrino irradiation, this behavior is not altered, for the amplitudes of neutrino luminosity 
of our interest. The irradiation tends to heat and push out the inf ailing matter. Hence, 
as the luminosity gets larger, the infall velocity gets smaller in the whole range and the 
solution curves go down. When the luminosity exceeds a critical value, there exists no point 
where the subsonic accretion flow joins with the free fall stream, via the Rankine-Hugoniot 
relations; the solution curves of the accretion flow Ur{r) never intersects with Mpost(^)- 

Prom the above explanation, we can also understand another important nature of the 
solutions; as the luminosity is raised with the mass accretion rate fixed, the radial derivative 
of the velocity of the accretion flow at the shock radius dur{r) / dr becomes larger for the inner 
solutions, whereas it becomes smaller for the outer solutions. If the flow is adiabatic, these 
derivatives coincide with that of dupos^{r)/dr at the shock radius for the critical luminosity. 
The neutrino heating, however, makes the former slightly smaller than the latter. This 
characteristics will be utihzed in section 4. 

Finally, let us discuss what happens when the luminosity exceeds the critical value. Since 
there is no steady solution, we should discuss this issue with time-dependent calculations. 
However, we can infer at least the initial response of the shock with a simpler discussion. 
We can find in figure 1 that when the luminosity exceeds the critical value, the downstream 
velocity at the shock wave becomes inevitably smaller than impost In addition, the down- 
stream entropy should be larger than that for the critical luminosity, because of the higher 
luminosity. We can regard this situation as a Riemann problem that has a discontinuity 
with a velocity slightly lower and an entropy slightly larger in the downstream than those 
satisfying the Rankine-Hugoniot relations. The solution of the Riemann problem tells us 
that, a shock wave will propagate outwards, while a weak rarefaction wave or another weak 
shock wave will propagate inwards at the same time (Courant and Friedrichs 1976); in our 
words, the shock wave revives. Thus, we can conclude that the condition for the shock to 
revive is that the neutrino luminosity exceeds the critical value. 
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3. Stability against Radial Perturbations 

In the previous section, we found that there are two branches of steady solutions for 
the accretion flow with a shock. In this section, we discuss which type of the solutions is 
realized, examining the stability of these solutions by global linear analysis. We consider the 
perturbations of neutrino luminosity, that is, non-adiabatic effects are taken into account. 

In this paper, only the stability against radial perturbations is considered. This is be- 
cause we did not take account of the effect of convection in obtaining the unperturbed flows 
in the previous section. There is, in general, a region where heating prevails over cooling 
under the irradiation of neutrinos; In this region, entropy increases in the direction of gravity 
and the flow is convectively unstable. Thus our solutions are convectively unstable against 
non-radial perturbations, which fact is of no interest to us here, though. 

Since we consider the perturbations accompanying variations of the positions of inner 
and outer boundaries, it is convenient to introduce a new variable x, defined as. 



x = ^^, (10) 
Ts - 

where denotes the shock radius. Then, the basic equations describing the time-dependent 
flows are written as. 



dp If dr^ .9r^1 dp 2pUr , p dur ^ Ur dp 

X— + [I - x)— > — -\ \ \ — = 0, (11) 



dt Ts — 1 dt dt \ dx r rg — Vi, dx Vg — dx 



dur 1 f ^ dry\ dur ^ Ur dur ^ 1 Idp ^ 

dt rs — r^X dt dt j dx rg — dx Vg — p dx 



de If dvg , , dr-y 1 (9e p dp 
X— + (1 - x)— > - — 



dt Ts — Ty \ dt dt J dx p^ dt 

dvg drApdp Ur de Ur p dp . 
^^ + (1-^)^ r^^ + ^ = (13) 



Tg — Ty \ dt dt J dx Tg — Ty dx Tg — Ty p^ dx 

Since we assume that the flow outside the shock is steady, the Rankine-Hugoniot rela- 
tions at the outer boundary are expressed as. 



p{Ur -Ug) ^ p{{Uf-Ug), 



(14) 
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p{Ur - Usf + p = Pf{U{ - UsY, (15) 

l,{ur - Usf + e + - = 1^{U{ - U^f , (16) 
2 P ^ 

where Us is the shock velocity, which can be written as 

In calculating the perturbations of heating and cooling rates, we assume, for simplicity, 
that the effective temperature of neutrinos is not perturbed. 

We perform a linear analysis, assuming that all the perturbed quantities have the time- 
dependence of exp(c(;t). In the rest of this section, the suffix 1 is attached to the perturbed 
quantities and the suffix to the unperturbed ones. Then, the equations describing the 
perturbations become 



d{l>i/p()) (){tlrl/ ll'ro) , rsO-TtyQ pi 2 
1 ^ \ LO \ 

OX UrO Po ro 



dx 
+ 



) 



Jso - _^ (9 In Po ^ dhiUro 



To 



^d\nUro 
ro ox 



dx 
1 91npo 



dx J \po 



Pi ^ Url 
UrO 



ro dx UrO dx 



Tsi ^ ,^ s r^i 

xvso ^ (1 - x)ruo — 

rso r^o 



0, (18) 



d{Url/Uro) ^ 1 



dx 
+ 

+ 



Ko Po 



+ 



Pmd{pi/po) , 1 4p Ro+pNo d{Ti/To) 

Po dx 
1 4pRo91nTopi 



u 



-UJ 



dx 



d\VLUro \ Url 



rO 



UrO 

1 pmd In Po 



UrO Po 



dx 



+ 



+ 



rO 

u \ din UrO 



dx J UrO UrO PO 9x 

1 16pRo + Pno ^ In To ^ Ti 

To 
2 



u 



Po dx 
2 1 pNo<91npo 



X 



UrO J dx 
XTsO \- (1 

rso 



roKo Po 



X Tyi 
x)r^o — 
r^o 



dx 



+ 



Po 



1 4pRo +PNo<91nTo 



To «rO 



Po 



dx 



GM 

2 2~ 
TqU^q 



TsO r^O 

TsO r„o 



0, 



(19) 
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Po dx po dx 

^ 4pRo+PNoQlnpo ^ 12pRo + (3/2)pno <91nTo | Uri 



Po dx Po dx j UrO 

+ i -[r^o - r^o) \ ^ ^ ^Po \ — 

Po UrO Po ox Po dx UrO Op ) po 

. 12pRo + (3/2)pno IGpRo+Pmdlnpo 
+ i ('^so -'^i.o) ^ — 

Po UrO Po OX 

_^ 48pRo + (3/2)pno ^InTp _ Tsq - r^o dq ^ 1 T\ 
Po dx UrO dT °/ To 



+ 



4pro+Pno CO dlnpo _ 12pRo + (3/2)pno dlnTp _ fjo^-r^dq 

Po UrO dx Po UrO dx UrO dr 

X <^ xr^o ^ (1 - x)ryo — 

I r^Q r„o 

q / r-si r^i\ r,o - r^p dq L^^i r,o - r^p dq ^ r^i _ .^f.. 

'^sO ^uO TT7 — J^i^eO~P 7^ — ^uO — U, [ZU) 

UrO \ fso r^o) ^ro dL^^ L^^o UrO dr^ r^o 
where, pro, Pno and rp are defined as 

PNo - (22) 

rp = XTso + (1 - a:)r^o. (23) 

All the variables, Uri/uro, Pi/po, Ti/To, rgi/rgp, r^i/r^o, L^^i/L^^o and an eigen-value uj can 
be complex. 

The outer boundary conditions {x = 1) become 



2 pfo Ufo UJ pfo \ Tsi /«rl , Pi , „ .r,A^ 
-2 ^ -TT^ PoOJTsO h PoMrO \ =0, (24) 

rtoPoUro r^^ Po J r^o \Uro Po, 



5 Pfo u}q 2u Pfo Ufo \ Tsi 

-?7^ 2~ 2a;rspppMr-o — 

2r 0^ Po ^^rO '^sO Po UrO J ^0 
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u T 
+2poM^o— + (Po^ro + Pm) — + (4pRo + Pm) ttt = 0, (25) 

UrO Po -to 



1 Mfo 1/2 MfO \ '^sl , 2 «rl , PRO Pi , 16pRo + 5/2pN0 „ ^OCN 

T 2- - ^^so t^rso-Uro ] — + u^o 4 + — = 0, 26) 

2rsO<o «rO / '^sO «ro Po Po Po To 

where pfo, Ufo are density and velocity of the free-fall at the shock in the unperturbed state. 

Since we set the neutrino-sphere as the inner boundary {x — 0) and the density is fixed 
there by definition, we impose the condition 

^ = 0, (27) 
Po 

at the inner boundary. In addition, we assume 

^ = 0, (28) 

UrO 

at the inner boundary, supposing that the infalling matter comes to rest on the neutron 
star. Since the neutrino temperature is assumed to be unperturbed, the variation of the 
luminosity is related with that of the neutrino-sphere radius as, 

Lu^ = 2^. (29) 

Above equations are solved as follows. We first take trial values of cu, r^i/ri,o, (rgi/rgo 
is set to be unity) and solve the conditions (24), (25), (26), and equations (18), (19), (20) 
from the outer boundary to the inner boundary. The solution obtained this way does not, 
in general, meet the inner boundary conditions (27), (28). Then we improve the value of a;, 
Tvi/fvo and repeat the procedure until (27), (28) are satisfied at the inner boundary. 

We found for a given unperturbed solution many modes with complex eigen-values 
(probably an infinite number of overtones, in fact) and only one mode with a real eigen-value. 
All modes with complex eigen-values arc damped, namely, the real part of the eigen-value is 
negative, owing to the thermal smearing effect of the neutrino irradiation. We show only the 
results for the mode with a real eigen-value. The growth rates are shown in figure 3. We can 
see that the inner solutions are always stable, while the outer solutions are always unstable. 
Further, when the luminosity reaches the critical value, both solutions become neutral. The 
eigen-functions are shown in figure 4. The ratios of the perturbation of luminosity to that of 
shock radius are shown in figure 5. Since the amplitude of luminosity perturbation is twice 
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as much as that of neutrino-sphere radius, we can see that the perturbations of inner and 
outer boundaries are in phase for the inner solutions, whereas they are in opposite phase for 
the outer solutions. 

We should note again that we investigated the stability only against radial perturbations. 
As several authors pointed out, spherical accretion shocks are, in general, unstable against 
non-radial perturbations (Houck & Chevalier 1992; Foglizzo 2002; Blondin et al. 2003). In 
this sense, we cannot say that the inner solutions are stable. Detailed studies of the stability 
against non-radial perturbations will be necessary with the effect of convection and neutrino 
irradiation properly taken into account. 



4. The Effects of Rotation 

In this section, we discuss the effects of rotation on the revival of shock. We extend 
the one dimensional analysis in section 2 to the two dimensional case. A brief summary of 
our results in this section has previously appeared (Yamada et al. 2004). Since as we found 
in the previous section, the radii of the outer solutions are unstable, we consider only the 
inner solutions in this section. The assumptions and simplifications mentioned in section 2 
are also adopted here except for the assumptions for the flows outside the shock, which will 
be mentioned later. Assuming axial symmetry, we use the spherical coordinates. Then the 
basic equations describing the steady flow inside the shock are written as, 

4|-(r>.) + -^^{sinepue) = 0, (30) 

dur , uedur uj + ul I dp GM 
or r Ou r p or r 

dug ugdue UrUe ulcotO 1 dp , , 

^^^ + ~^ + ^ ^' (32 

or r oB r r pr oB 

du^ ue du^ u^Ur upu^ cot B _ 
or r oB r r 

u (—- P_^\ ^^(^- p_dp\ ^ . 
^ \dr p^ dr J r \dB p^ dB J ' 

where u denotes velocity and the other notations are the same as those in section 2. 



Since we seek solutions with axial and equatorial symmetry, we impose the following 
conditions, 



dur dp dT 




(35) 



U(j, = U0 



de de de 



at 6* = 0, and 



dur du^ dp dT 




(36) 



de de de de 



at ^ = 7r/2. 



The outer and inner boundaries are set at the shock surface and neutrino-sphere, re- 
spectively. The neutrino-sphere will be oblate due to centrifugal force and, as a result, the 
neutrino flux will be anisotropic (Kotake et al. 2003). Such effect may have some influence 
on the explosion (Shimizu et al. 2001). Since our concern here is the hydrodynamical effects 
of rotation on the accretion flow, we assume that the neutrino-sphere is spherical and that 
the neutrino flux and temperature are isotropic. Unlike in section 2, the shock surface is not 
spherical and oblique to the flow, in general. 

As in the one dimensional calculations, we impose the condition that the density is 
10^^ g/cm^, at the inner boundary. The Rankine-Hugoniot relations for obhque shocks are 
imposed at the outer boundary. As for the flows outside the shock, we take the following 
approximations: i) the flows are radial except for rotation and free fall, ii) the density is 
independent of latitude, thus the accretion isotropic outside the shock, iii) the rotation 
frequency is independent of latitude and the specific angular momentum is conserved along 
each stream curve. Although these assumptions are artificial and not self-consistent, since 
the centrifugal force is small compared with the gravity or the inertia, the radial component 
of the velocity must be predominant in reality. 

We define the parameter / as the rotational frequency at the radius of 1000 km. We 
performed the calculations for / = 0.03 s~^ and 0.1 s~^, which correspond to the specific 
angular momenta averaged over the whole sphere outside the shock wave surface of 1.256 
and 4.187 x lO^^cm^/s, respectively. (The definition of the specific angular momentum is 
employed by Heger et al. (2000, 2003) in their calculations of stellar evolutions). 

In figure 6, we show the stream curves for the parameters M — 2.0 Mq/s, L^^ = 7.0-10^^ 
ergs/s, and / = 0.1 s^^. We can see that the fiow is pushed toward the equatorial plane 
owing to centrifugal force. The stream curves bend rather abruptly near the neutrino-sphere 
because the infalling matter is decelerated as it lands onto the neutron star surface and the 
radial velocity falls rapidly there. It is also noted that the shock surface is slightly deformed 
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into an oblate configuration, and the flows are directed toward the rotational axis at large 
radii. 

The angular dependence of the radial velocity inside the shock is shown in figure 7. At 
large radii, the velocity is larger near the rotation axis than the equatorial plane, whereas the 
opposite is true near the inner boundary. This behavior is attributed to the variation of cross 
section of the flow. The theory of Laval nozzle tells us that subsonic flows are decelerated 
when the cross section of nozzle gets larger, vice versa (Landau and Lifshitz 1983). As we 
saw in figure 6, the fiow is bent toward the equatorial plane near the ncutrino-sphcrc but 
toward the rotational axis near the shock. As a result, at small radii, the cross section of 
fiow becomes larger near the axis, whereas it gets smaller near the equatorial plane, and vice 
versa at large radii. Thus, we obtain the observed angular dependence of velocity. 

Next we discuss the effect of rotation on the critical luminosity. In figure 8, we show 
the dependence of the critical luminosity on the mass accretion rate. We can see that the 
smaller the mass accretion rate is, or the larger the rotation frequency is, the more reduced 
the critical luminosity is. According to the results of our calculations, if / is 0.03 s~^, and the 
accretion rate is 0.1 Mq/s, the critical value is about 25% smaller than that of the spherically 
symmetric fiow with the same mass accretion rate. If the rotation parameter is 0.1 s"^ the 
effect is more remarkable. For the mass accretion rates of 1.0 Mq/s and 0.5 Mq/s, we get 
the critical luminosities reduced by about 25 and 43% respectively, from the value for the 
spherically symmetric flow. 

Finally, we discuss the fate of the shock wave when the critical luminosity is reached. 
In section 2, we discussed the movement of the steady shock for the inner solutions; the 
radial derivative of the infall velocity increases when the neutrino luminosity is raised, and 
it becomes nearly equal to that of Upost{i^) when the luminosity reaches the critical value. 
The angular dependence of the velocity at the critical luminosity is shown in figure 7. We 
can see that the derivative is largest at = 0. This suggests that the steady fiow cease 
to exist first at the axis. Thus, we can infer that the stalled shock is revived and start to 
propagate outwards first at the rotational axis and that the jet-like explosion is followed, 
although time-dependent numerical simulations are necessary to confirm this. 

This break up of shock at the rotational axis and the resulting reduction of the crit- 
ical luminosity are also attributed to the changes of radial fiow velocities, which can be 
understood based on the theory of Laval nozzle, as we saw above. The variation of fiow- 
cross-section causes the acceleration of the radial fiow at large radii near the rotation axis 
whereas it leads to the deceleration near the equatorial plane. As a result, the revival of the 
stalled shock is more easily attained at the rotational axis for lower neutrino luminosities 
than in the spherically symmetric fiows. 
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5. Summciry and Discussion 

In this paper, we investigated the behavior of the stalled shock in the supernova core 
and the effects of rotation on its revival. The spherically symmetric accretion flows were re- 
investigated first. We found that there exist two branches of solutions when the luminosity 
of the irradiating neutrinos is below the critical value. The steady solution does not exist 
when the luminosity is higher than the critical value. 

In order to consider which solution is realized below the critical luminosity, we then 
examined the stability of the spherical accretion flows against radial perturbations, and 
found that the inner solutions are stable while outer solutions are unstable. In fact, for 
small luminosities expected for the early phase of the shock stagnation, the radii of the 
outer solutions are extremely large (flgure 2). Thus the outer solutions will not be realized 
in reality. The inner solutions are always stable until the neutrino luminosity reaches the 
critical value, where the two branches merge and the solution becomes neutrally stable. 
It should be noted, however, that the stability of the outer solutions may be affected by 
rotation, which will be a future work. 

We did not discuss the stability against non-radial perturbations, since our model are 
convectively unstable. Several authors showed the spherical accretion shock is unstable 
against non-radial perturbations. Blondin et al. (2003) demonstrated the instability for 
adiabatic perturbations by numerical simulations. Houck & Chevalier (1992) examined such 
a stability taking into account the cooling processes whose rates are determined locally. In 
the present situation, however, the heating rates are not determined locally because the 
flow is irradiated by neutrinos coming out of the proto-neutron star; the perturbations of 
the flow affect the neutrino luminosity, and the latter then affects the former, by varying 
the heating and cooling rates. Thus, the problem is different from those they treated. The 
detailed analysis in such a situation is a future work. Furthermore, the stability analysis 
with rotation taken into account is also an interesting issue. 

Finally, we discussed the effects of rotation on the revival of the stalled shock. We 
showed that rotation lowers the critical luminosity. However, we claimed that the shock 
revival would first take place at the rotation axis. We took the rotation frequency of 0.03 
and 0.1 s~^ at the radius of 1000 km and found that the rotation effect is substantial. The 
rotation frequency of newly born neutron stars guessed from the spin down rate of young 
pulsars is about several ten per second (Kaspi et al. 1994; van der Swaluw and Wu 2001). 
Thus the frequency we adopted here is one order of magnitude higher than the value inferred 
from the observations of young pulsars. It might be possible, however, that this discrepancy 
is removed by some mechanisms such as magnetic breaking in the late phase of explosion, 
which could take away the excessive angular momentum from the proto-neutron stars. In 
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fact, the values we adopted are consistent with those obtained in calculations of the evolution 
of rotating stars (Heger et al. 2003). Anyway, the angular momentum distribution inside 
massive stars are not yet well understood because of our poor knowledge on the mechanism 
of angular momentum transform (cf. Heger et al. (2000, 2003)). If the angular momentum 
of the infalling matter is large enough, our results suggest that rotation assists the shock 
revival. 

In this paper, we employed various assumptions and simplifications in the formulations. 
We assumed that the neutrino-sphere is spherical and the neutrino fiux is isotropic. If the 
collapsing star is rotating, the proto-ncutron star will have angular momentum. Kotakc ct 
al. (2003) investigated how the neutrino-sphere is affected by rotation. They showed that 
the neutron star becomes oblate and, as a result, the emerging flux is more concentrated 
to the rotation axis. One the other hand, Shimizu et al. (2001) studied how the explosion 
is affected when the neutrino flux is greater near the rotation axis, and claimed that such 
neutrino anisotropy can trigger a jet-like explosion in otherwise failed models. Such effects 
were neglected in this paper just for simplicity and are expected to further reduce the critical 
luminosity, which we will study in the future. 

The convection should be also affected by rotation. Fryer & Heger (2000) studied such 
effects by numerical simulations. They found that the centrifugal force reduce the effective 
gravity, and, thus, the efficiency of convection is lowered near the equatorial plane. As a 
result, the explosion energy was also decreased. The explosion, however, becomes jct-fikc. In 
order to investigate the effect of convection in the time-independent models, it is necessary 
to introduce some simplification, such as a mixing length theory or other formulation to 
smooth the entropy gradient. This is currently being undertaken. 

Finally the collapsing stars may have magnetic fields, which are then amplified by ro- 
tation. The recent works suggest that the magnetic fields arc amplified by the differential 
rotation or magneto-rotational instability (Akiyama et al. 2003; Kotake et al. 2004; Ardeljan 
et al. 2004) . If the magnetic energy is amphfled to be comparable to the rotational energy, it 
must affect the accretion flow and the behavior of the stalled shock as they suggested. This 
issue is also a future work. 
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Fig. 1. — a. Inner solution curves for spherical steady accretion flows with M = 2.0 Mq/s, 
Lj/g = 5, 6, 7, 8 • 10^^ ergs/s and 8.3167- 10^^ ergs/s (critical value) from left to right, b. Outer 
solution curves for steady accretion flows with M = 0.2 Mq/s, L,y^ = 6, 7, 8 • 10^^ ergs/s and 
8.3167 • 10^^ ergs/s from left to right. Dashed hues denote shock jumps. Dotted lines show 
the downstream values satisfying the Rankine-Hugoniot relations at each radius. 
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Fig. 2. — Radii of shock surface. Solid curves denote the inner solutions for M = 
1.0,2.0,3.0,4.0 Mq/s from left to right. Dashed curves display the outer solutions for 
M = 1.0,2.0,3.0,4.0 Mq/s from left to right. The crosses represent the critical points. 
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Fig. 3. — Growth rates for radial perturbations. Solid curves show the inner solutions for 
M = 1.0,2.0,3.0,4.0 Mq/s from left to right. Dashed curves depict the outer solutions for 
M = 1.0, 2.0, 3.0, 4.0 Mo/s from left to right. 
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Fig. 4. — a. eigen-functions of radial perturbations for the inner solutions with M = 2.0 
Mq/s, L^^ = 6, 7, 8 • 10^^ ergs/s and 8.3167 • 10^^ ergs/s (critical value) from left to right, 
b. eigen-functions of radial perturbations for the outer solutions with M = 0.2 Mq/s, 
L^^ = 6, 7, 8 ■ 10^2 ergs/s and 8.3167 ■ 10^^ gj-gg/g fj-^j^ ^ight to left near the shock. Note that 
the functions are defined in different ranges of radius for different models. 
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Fig. 5. — Ratio of perturbation of luminosity to that of shock radius. Sohd curves denote 
the inner solutions for M = 1.0, 2.0, 3.0, 4.0 Mq/s from left to right. Dashed curves present 
the outer solutions for M = 1.0, 2.0, 3.0, 4.0 Mq/s from left to right. 




Fig. 6. — Meridian section of a collapsing star and stream curves for M = 2.0 Mq/s, L^,^ — 
7 ■ 10^^ ergs/s and rotation frequency 0.1 s~^ at r = 1000 km. 
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Fig. 7. — Radial velocities at latitude 9 = 0, 7r/8, 7r/4, 37r/8 and n/2 for the rotating steady 
accretion flow with M = 2.0 Mq/s, = 8.10 • 10^^ ergs/s (critical value) and rotation 
frequency 0.1 s^^ at r = 1000 km from top to bottom. Dotted line shows the downstream 
value satisfying the Rankine-Hugoniot relations at each radius. 
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Fig. 8. — Critical luminosity. Solid curve denotes the spherically symmetric case, and the 
others correspond to the rotational models. The rotation frequencies are 0.03, 0.1 s~^ at 1000 
km for dotted and dashed curves, respectively. No steady solution exists for the luminosity 
larger than the critical value. 



